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Ordering of the Heisenberg spin glass with the nearest-neighbor Gaussian coupling is investigated 
ff) ■ by equilibrium Monte Carlo simulations in four and five dimensions. Ordering of the mean-field 

Heisenberg spin glass is also studied for comparison. Particular attention is paid to the nature 
of the spin-glass and chiral-glass orderings. Our numerical data suggest that, in five dimensions, 
£Nj . the model exhibits a single spin-glass transition at a finite temperature, where the spin-glass order 

accompanying the simultaneous chiral-glass order sets in. In four dimensions, the model exhibits a 
marginal behavior. Chiral-glass transition at a finite temperature not accompanying the standard 
spin-glass order is likely to occur, while the critical region associated with the chiral-glass transition 
is very narrow suggesting that the dimension four is close to the marginal dimensionality. 
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INTRODUCTION 



m . 

: 

In numerical studies of spin glasses (SGs), much effort has been devoted to clarify the properties of the so-called 
Edwards- Anderson (EA) modeli. Most of these numerical works on the EA model have concentrated on the Ising EA 
model. It is also very important, however, to clarify the properties of the corresponding Heisenberg model. This is 
i t ' simply due to the fact that many of real SG magnets are Heisenberg-like rather than Ising-like in the sense that the 
magnetic anisotropy is considerably weaker than the isotropic exchange interaction^** 

Indeed, several numerical works have been performed on the Heisenberg EA model. Earlier numerical studies 
I . suggested that, in apparent contrast to experiments, the isotropic Heisenberg SG in three dimensions (3D) did not 
exhibit an equilibrium SG transition at any finite temperature. This observation leads to general belief 

that the weak random magnetic anisotropy is crucially important in realizing a finite-temperature SG transition and 
q ■ a stable SG phase, which causes a crossover from the T g = isotropic Heisenberg behavior to the T g > anisotropic 
'— Ising behavior. The expected Heisenberg-to-Ising crossover, however, has not been observed experimentally, and this 
■ puzzle has remained unexplained^^ 

Meanwhile, a novel possibility was suggested by one of the present authors (H.K.) that the 3D Heisenberg SG might 
. exhibit an equilibrium phase transition at a finite temperature, not in the spin sector as usually envisaged, but in the 
chirality sector, i.e., might exhibit a chiral-glass transition. 6 Chirality is a multispin variable representing the sense 
or the handedness of local noncoplanar spin structures induced by spin frustration. In the chiral-glass ordered state, 
the chirality is ordered in a spatially random manner while the Heisenberg spin remains paramagnetic. References 
&2i2ii2iiiii2 claimed that the standard SG order associated with the freezing of the Heisenberg spin occurred at a 
i^: 1 temperature lower than the chiral-glass transition temperature at T — Tsg < Tqg, quite possibly Tsg = 0. It means 
that the spin and the chirality are decoupled on long length scales (spin-chirality decoupling). In fact, based on 
such a spin-chirality decoupling picture, a chirality scenario of the SG transition has been advanced, which explains 
the experimentally observed SG transition as essentially chirality driven 6 ^. Note that the numerical observation of 
""j" 1 , a finite-temperature chiral-glass transition in the 3D Heisenberg SG of Refsi&SiSiiSiiiii^ is not inconsistent with the 
' earlier observations of the absence of the conventional SG order at any finite temperature. 

Recently, however, in a series of numerical studies on the 3D Heisenberg EA model, Tohoku group criticized the 
earlier numerical works, claiming that in the 3D Heisenberg SG the spin ordered at a finite temperature and that the 
SG transition temperature might coincide with the chiral-glass transition temperature, i.e., Tsg — ^CG > 0^2*14. By 
^ • contrast, Hukushima and Kawamura maintained that in 3D the spin and the chirality were decoupled on sufficiently 
k^j ' long length scales, and that Tsg < Tea—, supporting the earlier numerical results. The situation in 3D thus remains 
. controversial. 

Under such circumstances, in order to shed further light on the nature of the ordering in 3D, it might be useful 
to study the problem for the general space dimensionality D, particularly for dimensions higher than D = 3. In the 
limit of infinite dimensions D — > oo, the model reduces to the corresponding mean- field model, i.e., the Heisenberg 
Sherrington-Kirkpatrick (SK) model. In the case of equal weights of the ferromagnetic and antiferromagnetic inter- 
actions, the SK model is known to exhibit a single continuous SG transition. Hence, in the D — > oo limit, the order 
parameter of the transition is the Heisenberg spin itself, with no exotic phase such as the chiral-glass phase. Further- 
more, the SG ordered state of the SK model is known to exhibit a hierarchical type of replica-symmetry breaking 
(RSB), i.e., a full RSB. 

Then, questions which naturally arise are: (i) What is the lower critical dimension (LCD) of the SG order cq G ? (ii) 
Is df G the same as the LCD associated with the chiral-glass order d c ? 
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Concerning the point (i), several earlier numerical studies including the high-temperature expansion^ and the 
numerical domain-wall renormalization-group calculation ££L suggested that <if G might be close to four. Meanwhile, 
Anderson and Pond argued that df G = First Monte Carlo (MC) simulation on the high-dimensional Heisenberg 
EA model was performed by Stauffer and Binder—. By studying the temporal decay of the EA order parameter, 
they suggested that a finite-temperature SG order occurred in D = 5 and 6, but not in D < 4 19 . More recently, 
the 4D Heisenberg EA model was studied by Coluzzi by equilibrium MC simulation 20 . By examining the behavior 
of the Binder ratio, she suggested the occurrence of a finite-temperature SG transition in D = 4 in contrast to the 
suggestion of Ref— . There seems to be no consensus as to the point (i) . 

The point (ii) above is closely related to the controversy regarding whether the spin-glass and the chiral-glass 
orders occur simultaneously or separately in 3D. To the authors' knowledge, concerning the chiral-glass order in 
D > 4 dimensions, no calculation has been reported so far. In the present paper, we wish to fill this gap. We study 
both the spin-glass and the chiral-glass orders of the Heisenberg EA model in both 4D and 5D by means of a large- 
scale equilibrium MC simulation. In particular, we simulate larger lattices and lower temperatures than those covered 
in RefiS .. For comparison, a simulation is also performed on the mean-field Heisenberg SK model corresponding to 
D = oo. 

Our data suggest that, in 5D, the model exhibits a single SG transition at a finite temperature reminiscent to the 
one of the Heisenberg SK model. The chirality orders simultaneously with the spin, but it behaves as the composite 
operator of the spin, not as the order parameter. The SG ordered state in 5D accompanies a peculiar type of RSB, 
most probably a one-step-likc RSB, which is different in character from the full RSB realized in the SK model. In 4D, 
by contrast, the model exhibits a pure chiral-glass transition at a finite temperature, not accompanying the standard 
SG order. The critical region associated with the chiral-glass transition, however, is very narrow, suggesting that the 
4D model lies close to the marginal dimensionality. The chiral-glass ordered state accompanies a one-step-like RSB. 

The present paper is organized as follows. In iflTl we introduce our model and explain some of the details of the 
MC calculation. Various physical quantities calculated in our MC simulation are defined in £11111 The results of our 
MC simulation on the 4D, 5D and SK models are presented in fllVI SectionElis devoted to summary and discussion. 

II. THE MODEL AND THE METHOD 

The model we consider is the isotropic classical Heisenberg model on a 4D or 5D hypercubic lattice, with the 
nearest-neighbor Gaussian coupling. The Hamiltonian is given by 

H = - 2J JijSi ■ Sj , (1) 

<ij> 

where Si = (Sf , Sf , Sf) is a three-component unit vector, and < ij > sum is taken over nearest-neighbor pairs on 
the lattice. The nearest-neighbor coupling is assumed to obey the Gaussian distribution with a zero mean and a 
variance J 2 . 

For comparison, we also simulate the corresponding infinite-ranged model, i.e., the Heisenberg SK model corre- 
sponding to D — ► oo. In the SK model, the Gaussian coupling J^- works between all possible pairs of total iV spins 
with a zero mean and a variance J 2 /N . 

We perform an equilibrium MC simulation on these models. In 4D, the lattices studied are the hypercubic lattices 
with N = L 4 sites with L = 4, 6, 8, and 10, whereas in 5D, N = L 5 with L = 3, 4, 5, 6 and 7. In the case of the SK 
model, N is taken to be N = 32, 64, 128, 256 and 512. In all cases, we impose periodic boundary conditions in all D 
directions. The sample average is taken over 64-200 independent bond realizations, depending on the system size L 
and the lattice dimensionality D. Error bars of physical quantities are estimated by the sample-to-sample statistical 
fluctuation over the bond realizations. 

In order to facilitate efficient thermalization, we combine the standard heat-bath method with the temperature- 
exchange technique—. Care is taken to be sure that the system is fully equilibrated. Equilibration is checked by the 
following procedures. First, we monitor the system to travel back and forth many times along the temperature axis 
during the the temperature-exchange process (typically more than 10 times) between the maximum and minimum 
temperature points. We check at the same time that the relaxation due to the standard heat-bath updating is 
reasonably fast at the highest temperature, whose relaxation time is of order 10 2 Monte Carlo steps per spin (MCS). 
This guarantees that different parts of the phase space are sampled in each "cycle" of the temperature-exchange run. 
Second, we check the stability of the results against at least three times longer runs for a subset of samples. Third, we 
use the method recently developed in Refsi 22 ' 23 for the Gaussian coupling, in which a certain quantity is calculated in 
two ways, each of which is expected to approach the asymptotic equilibrium value either from above or from below. 
Further details of our MC simulation are given in Table [I] 
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TABLE I: Detailed conditions of the MC simulation. Here, D represents the spatial dimensionality, N the total number of 
spins, -/V samp the total number of samples, Nt the total number of temperature points used in the temperature-exchange run, 
Imax/J and T m i n /J the maximum and minimum temperatures in the temperature-exchange run. 



D 


N 




N T 


^max/ J 


Txaixi/ J 




4 4 


140 


20 


1.0 


0.16 


4 


6 4 


100 


34 


0.65 


0.21 




8 4 


80 


50 


0.58 


0.21 




10 4 


80 


50 


0.48 


0.256 




3 5 


120 


24 


1.2 


0.24 




4 5 


100 


36 


1.0 


0.3 


5 


5 5 


96 


46 


0.8 


0.31 




6 5 


96 


56 


0.75 


0.40 




7 5 


64 


50 


0.75 


0.537 




32 


200 


30 


0.5 


0.1 




64 


200 


30 


0.5 


0.1 


oo(SK) 


128 


200 


30 


0.5 


0.1 




256 


100 


34 


0.5 


0.1 




512 


100 


34 


0.5 


0.1 



III. PHYSICAL QUANTITIES 



In this section, we define various physical quantities calculated in our simulation below. 



A. Spin-related quantities 



By considering two independent systems ("replicas") described by the same Hamiltonian one can define an 
overlap variable. The overlap of the Heisenberg spin is denned as a tensor variable between the fi and v components 
(fi, u=x,y,z) of the Heisenberg spin, 

1 N 

q ^ = N S il )S i» ' V = X ' ^' Z ) ' ( 2 ) 

i=l 

where and are the i-th Heisenberg spins of the replicas 1 and 2, respectively. In our simulation, we prepare the 
two replicas 1 and 2 by running two independent sequences of systems in parallel with different spin initial conditions 
and different sequences of random numbers. In terms of these tensor overlaps, the SG order parameter is defined by 

qP = [tf s )} , q i = C . ( 3 ) 



fi,v—x,y,z 



while the associated spin Binder ratio is defined by 



where (• • •) represents the thermal average and [• • •] the average over the bond disorder. Note that g s is normalized here 
so that, in the thermodynamic limit, they vanish in the high-temperature phase and gives unity in the nondegenrate 
ordered state. 

The spin-overlap distribution function is originally defined in the tensor space, since the relevant spin-overlap has 
3x3 = 9 independent components. For the convenience of illustration, we introduce here the diagonal spin-overlap 
distribution function, 

-Ps^diag) = [(<%diag ~ Sdiag))] , (5) 
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defined in terms of the diagonal overlap Qdiag which is the trace of the tensor overlap gw's, 

1 N 

<7diag = q ^ = ' ■ ( 6 ) 

Note that the spin Binder ratio Eq.Q is defined from the full tensor overlap distribution function, but cannot be 
derived solely from the diagonal overlap distribution function Eq.©. 

In zero field, the distribution function -P s (<Zdiag) is symmetric with respect to Qdiag = 0. In the high-temperature 
phase, each q^ v (p, v = x, y, z) is expected to be Gaussian-distributed around q^ v — in the L — > oo limit, and so is 
qdiag- In the hypothetical SG ordered state, reflecting the fact that gdiag transforms nontrivially under the independent 
global 0(3) spin rotations on the two replicas, even a self-overlap part of P s (gdiag) develops a nontrivial shape, i,e„ it 
is not a simple delta-function located at the EA SG order parameter ±qf A . Hence, we need to clarify first how the 
function P s (gdiag) looks like in the possible SG ordered state with a nonzero ±qf A — . 

Let us hypothesize here that there exists a spm-glass ordered state characterized by a nonzero EA SG order 
parameter qf A > 0. One can show in the L — > oo limit that the self-overlap part of -P s (<Zdiag) is given by 

Pin \- 3 ^ I I ~ gdiag / gP + gdiag \ ( „s 

which is illustrated in FigQ] The derivation of Eq.Q is given in the Appendix. Note that diverging (5-function peaks 
appear at <?diag = ±3"<? EA , n °t at Qdiag = ±<? EA - If the SG ordered state accompanies RSB, the associated nontrivial 
contribution would be added to the one given by Eq.Q. In any case, an important observation here is that, as long 
as the ordered state possesses a finite SG long-range order (LRO), the diverging peak should arise in -P s (9diag) a t 

L 1 „EA 

<?diag - ±3<? s ■ 




FIG. 1: Sketch of the self-overlap part of the diagonal spin-overlap distribution function P s (gdia g ) in the thermodynamic limit, 
expected when there exists a finite SG long-range order with a nonzero qf A > 0. 



B. Chirality-related quantities 



We define the local chirality at the i-th site and in the /i-th direction, Xi^j for three neighboring Heisenberg spins 
by the scalar 

Xip. = Si+e„ ■ (Si X 5,;_e M ) , (8) 

where e M (fx = Xi, x%, ■ ■ ■ , xd) denotes a unit vector along the /i-th axis. By this definition, there are in total DN 
local chiral variables. 
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In the case of the SK model, since there is no "neighbors" in the model, the definition of the local chirality 
accompanies some difficulties. Here, just for convenience, we number the N spins arbitrarily, and define the local 
chirality by 



Xi = S i+1 ■ { S i x Si-l) ■ 

Then, there are in total N chiral variables. 

The mean local amplitude of the chirality, x, may be defined by 



(9) 



X = 



1 N 

[(Xi M )] ' 



(10) 



= 1 ll=Xl,X2,---,XD 



for D = 4 and 5, and by 



\ 



1 N 



(ii) 



for the SK model. Note that the magnitude of x tells us the extent of the noncoplanarity of the local spin structures. 
In particular, this quantity vanishes for any coplanar spin configuration. 

As in the case of the Heisenberg spin, one can define an overlap of the chiral variable by considering the two replicas. 
In the cases of D = 4 and 5, it is defined by 



1 N 



(1) (2) 



(12) 



i=l ^=xi,X2,---.x D 



where Xu? an d Xu? represent the chiral variables of the replicas 1 and 2, respectively. In the case of the SK model, 
it is defined by 



1 N 

J- (!) (2) 



N 



i=i 



In terms of this chiral overlap q x , the chiral-glass order parameter is defined by 
The associated chiral-glass susceptibility may be defined by 



Xx = DN[(q 2 x )] , 

in the cases of D = 4 and 5, while in the case of the SK model, it is defined by 

X x = N[{q x )\ . 



(13) 

(14) 
(15) 
(16) 



Unlike the spin variable, the local magnitude of the chirality is somewhat temperature dependent. In order to take 
account of this effect, we also consider the reduced chiral-glass susceptibility \x by dividing \x by the appropriate 
powers of x, 



x x -- 4 



The Binder ratio of the chirality is defined by 



9x = o 3 



The distribution function of the chiral overlap q x is defined by 

P x {q' x ) = Mq' x -q x ))] . 



(17) 



(18) 



(19) 
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IV. MONTE CARLO RESULTS 

In this section, we present our MC results on the Heisenberg EA models in 4D, 5D, and in the SK limit (corre- 
sponding to D = oo). 

A. Chiral Binder ratio 

In Fig|2 we show the Binder ratio of the chirality in the cases of (a) 4D, (b) 5D and (c) the SK model, respectively. 
In all these cases, g x exhibits a negative dip, while its temperature and size dependence is somewhat different from 
each other. 

In 4D, with increasing the lattice size L, the negative dip tends to deepen while the dip temperature T^- lp is almost 
kept constant at around T/J — 0.38: See the inset of FigEI One can also see from FigElthat g x for various L cross 
in the negative region of g x at temperatures slightly above Tdi p . The occurrence of a negative dip deepening with L, 
accompanied by a crossing on the negative side of g x , is similar to the one previously observed in the corresponding 
3D model 1011 ' 12 , although, in 3D, Tdip tends to shift toward lower temperature with increasing L. As argued in 
Ref4Sill*l£, the occurrence of a negative dip deepening with L is a signature of the occurrence of a phase transition 
in the chiral sector. By making a linear extrapolation of T^i p (L) with respect to as shown in the inset of Fig[21 
(a), we estimate the bulk chiral-glass transition temperature as Tqg/J = 0.38(2). Below Tec the curves for L > 6 
almost merge into a curve exhibiting the nontrivial temperature dependence. Such a behavior suggests that the chiral 
ordered state accompanies a nontrivial phase-space structure, i.e., RSB. 

In 5D, although the negative dip tends to deepen up to the size L = 5, it tends to become shallower for L > 6. In 
contrast to the 3D and 4D cases, the dip temperature Tdip tends to shift toward higher temperature with increasing L. 
The observed temperature and size dependence of g x strongly suggests again that the limit Tdi P (L — > oo) corresponds 
to a transition temperature in the chiral sector. The 1/L-extrapolation of Tdi P (L) to L — > oo, shown in the inset of 
Fig0 (b), yields the estimate T CG /J = 0.62(2). 

In the SK case, the negative dip of g x becomes shallower with increasing L, and Tdi p shifts toward higher temper- 
ature. In the inset, T dip (N) is plotted as a function of jV -1 /325i2£, wn j cn y i e lds T CG /J = 0.31(2). The estimated 
chiral-glass ordering temperature agrees within errors with the exactly-known SG transition temperature of the Heisen- 
berg SK model, T$q/J = 1/3. This coincidence simply confirms the fact that, at the SG transition of the Heisenberg 
SK model, Heisenberg spins order in a noncoplanar manner, which necessarily accompanies the onset of a nonzero 
chiral-glass LRO. As is evident, in the SK case, the order parameter of the transition is not the chirality, but the 
Heisenberg spin itself. 

B. Chiral autocorrelation function 

More direct measure of the chiral-glass transition may be obtained from the equilibrium dynamics of the model. 
We compute the autocorrelation function of the chirality defined by 

1 N 

C X W = M E E [(XiMxi^t + to))] , (20) 

i— 1 ^—Xi,X2,---,Xe> 

where the "time" t is measured here in units of MCS. In computing l|20() . the simulation is performed according to 
the standard heat-bath updating without the temperature-exchange procedure, while the starting spin configuration 
at t = to is taken from the equilibrium spin configurations generated in our temperature-exchange MC runs. 

We show in Fig|3]the time dependence of C x {t) on a log-log plot for the cases of (a) 4D, and (b) 5D. To check the 
possible size dependence, the data for the two largest lattice sizes are given together, one denoted by symbols and the 
other by thin lines. In the chiral-glass ordered state with a nonzero q x A , C x {t) in the L — + oo limit should exhibit an 
upward curvature, tending to q^ A > 0. In the disordered phase, C x (t) should exhibit a downward curvature. Just at 
T = Tqq, the linear behavior corresponding to the power-law decay is expected. As shown in the figures, in the time 
region where the finite-size effect is negligible, C x (t) shows either a downward curvature characteristic of the disordered 
phase, or an upward curvature characteristic of the ordered phase, depending on whether the temperature is higher or 
lower than a critical value. In 4D, the chiral-glass transition temperature estimated in this way is Tqq/ J — 0.38(2), 
while in 5D it is Tqg/J = 0.60(2). Both are close to our estimate above based on the chiral Binder ratio. Our 
observation that C x (t) exhibits an upward curvature below Tqg indicates that the chiral-glass ordered state is "rigid" 
with a nonzero long-range order parameter. 
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From the above analysis, the occurrence of the chiral-glass LRO in 4D and 5D seems now well established. The 
next question is whether the chiral-glass order accompanies the standard SG order. If the SG order occurs at the 
same temperature as the chiral-glass order, the transition is likely to be of the standard type, at least in the sense 
that the order parameter of the transition is the spin, not the chirality. (Here, recall that the SG ordered state in the 
Heisenberg SG inevitably accompanies the chiral-glass order as long as the spin is frozen in a noncoplanar manner.) 
By contrast, if the SG order occurs at a temperature below the ordering temperature of the chirality, it means the 
unusual situation, i.e., the occurrence of the spin-chirality decoupling, the pure chiral-glass transition and the pure 
chiral-glass ordered state. To clarify this issue, we examine the spin Binder ratio in the next subsection. 



C. Spin Binder ratio 

In FigQ] we show the temperature and size dependence of the Binder ratio of the spin in the cases of (a) 4D, (b) 5D 
and (c) the SK model. In each figure, magnified figure is embedded to show the detailed behavior of g s in the region 
of interest. The arrow in each figure indicates the location of the transition point of the chirality determined above. 
In 4D, as can be seen from Fig0](a), g s for the range of sizes 4 < L < 8 appears to almost merge at a temperature 
T/J ~ 0.4 close to the chiral-glass transition temperature determined above. This seems to suggest that the spin 
sector also becomes critical at T ~ Tcg> which may indicate the simultaneous spin-glass and chiral-glass transition 
at T/J ~ 0.4. Quite remarkably, however, the spin Binder ratio g s for our largest size L = 10 comes definitely below 
the curves for L < 8. In fact, the L — 10 curve lies below the L < 8 curves more than four sigmas, and the observed 
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FIG. 3: Temporal decay of the equilibrium chiral autocorrelation function at various temperatures in 4D (a), and in 5D (b). In 
(a), symbols represent the data of L = 8 and thin lines those of L = 10. In (b), symbols represent the data of L — 6 and thin 
lines those of L — 7. Solid straight lines represent the power- law fits of the data at T/J — 0.38 (a), and at T/J = 0.60 (b). 



departure is statistically well significant. From this observation, we conclude that, although the spin sector once 
becomes near-critical at the chiral-glass transition T = Tcg 011 short length scales, it eventually remains off-critical 
(paramagnetic) at T — Tqg on longer length scales. Namely, the spin-chirality decoupling previously observed in 3D 
in Refpi2> seems to come into play in 4D as well, which keeps the spin sector being paramagnetic even below T = Tqg- 
If this is the case, the transition in 4D is a pure chiral-glass transition and the ordered state is a pure chiral-glass state 
not accompanying the standard SG order. Even in certain temperature range below Tcg, <?s is expected to approach 
zero in the L — > oo limit, which seems consistent with the present data. We interpret the strange structure observed 
in g s near T = Tcg for larger L as a remanence of the near-critical behavior of the spin at the chiral-glass transition. 
Meanwhile, due to the lack of our data in the lower temperature region T / J 5 0.2, it is difficult to determine from the 
present data whether the SG transition occurs either only at zero temperature, Tsg = 0, or at a finite temperature 
below the chiral-glass transition temperature, < Tsg < Tcg- Nevertheless, FigQJa) strongly suggests that Tsg, if 
it is nonzero, is less than 0.2 J. 

We note in passing that the spin Binder ratio g s of the same 4D model was calculated by Coluzzi for smaller 
sizes L = 3,4 and 5 and at higher temperatures T/J > 0.521. She observed that g B for L — 3,4 and 5 appeared to 
merge around T/J ~ 0.5, and suggested that there occurred a standard SG transition at T/J ~ 0.5. In the present 
calculation made for larger lattice sizes and for lower temperatures, although we indeed observed a near-merging 
behavior of g s , it occurred at a temperature T/J ~ 0.4 near Tcg, somewhat lower than the estimate of Refi**, and 
most importantly, g s for larger lattices show clear deviation from the merging behavior, suggesting that the Heisenberg 
spin remains paramagnetic even below the chiral-glass transition point. 

Now, we turn to the spin Binder ratio in 5D. As shown in FigQ](b), g s for 3 < L < 7 show a crossing at a temperature 
around T/J ~ 0.60, strongly suggesting that the SG order occurs at Tsg/ J — 0.60(2). At lower temperatures, g s for 
larger L tends to come down again, exhibiting a behavior reminiscent to the one observed in 4D. In particular, at low 
enough temperatures, g s decreases with L. The observed non-monotonic temperature dependence and the peculiar 
size dependence of g s below Tsg suggests that the SG state below Tsg might be a nontrivial one accompanied by a 
peculiar RSB, possibly the one-step-like one as observed in the chirality sector. 

In the SK case, the calculated g s exhibits a clear crossing behavior at around T/J = 1/3, an exactly known SG 
transition temperature of the model. At lower temperatures, g s monotonically increases with N in contrast to the 4D 
and 5D cases, and eventually appear to converge to the temperature-dependent values less than unity. The observation 
that the asymptotic g s (L — > oo) in the region < T < ~J = Tsg differs from unity reflects the fact that the SG 
ordered state of the SK model accompanies a full (or hierarchical) RSB. 
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FIG. 4: Temperature and size dependence of the spin Binder ratio g s in 4D (a), in 5D (b), and in the SK model (c). Magnified 
figures are given as insets. The arrow in each figure represents the location of the chiral-glass transition point. 



D. Chiral overlap distribution function 



In FigEl we show the overlap distribution function of the chirality, P x (q x ), well below Tqg, at around T ~ §7cg- 
In 4D, in addition to the "side peaks" corresponding to q x = ±q x A which grow and sharpen with increasing L, a 
"central peak" appears at q x = for L > 6, which also sharpens and gets higher with increasing L. These features 
were reminiscent to the ones observed in 3D, which was interpreted as a signature of the one-step-like RSB in the 
chiral-glass ordered statei2*ii*i£. Our present data suggest that the chiral-glass state in 4D also accompanies the 
one-step-like RSB as in the 3D case. This is fully consistent with the observed behavior of the chiral Binder ratio. 

A central peak in P x (q x ) is also observed in 5D in the largest lattice size L — 6. This suggests the occurrence of 
the one-step-like RSB also in 5D. 

In the SK case, the calculated P x {q x ) exhibits the side peaks at q x = ±q x A only, without a central peak at q x = 
for any size studied. Instead, the value of P x (0) gradually decreases with increasing N, where 1/iV-extrapolation of 
P x (0) to N — > oo gives a nonzero value, P x {q x = 0, N = oo) ~ 6.0 x 1(T 4 . Thus, in the SK case, the chirality exhibits 
the standard full RSB, in apparent contrast to the 4D and 5D cases. 



E. Spin overlap distribution function 



In FigHjl we show the diagonal spin-overlap distribution function, Ps(<?diag) 5 for the cases of (a) 4D, (b) 5D and (c) 
the SK model. The temperatures are taken to be the same as those for the chiral overlap distribution shown in FigO 
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FIG. 5: Overlap distribution function of the chirality, in 4D at T/J = 0.26 (a), in 5D at T/J = 0.40 (b), and in the SK model 
at T/J = 0.10 (c). 



i.e., about §Tcg- 

In the SK case shown in FigEJc), the shape of -Ps('Zdiag) is similar to the one of Fig^ with symmetric diverging 
peaks observed at <?diag — ±0.2. These peaks are then identified with the ones expected at ±^qf A when there is a 
finite SG LRO. This observation is fully consistent with the fact that the standard SG LRO with a nonzero qf A is 
realized in the ordered state of the SK model. In 5D, the growing symmetric peaks also appear at (faiag — 

±0.18 for 

L > 5, suggesting that the SG LRO characterized by a nonzero qf A is realized. 

In 4D, by contrast, no peaks corresponding to ±^qf A are observed, at least within the range of sizes we simulate. 
Instead, P s (<?diag) exhibits a marginal behavior, staying nearly flat with a plateau-like structure at |(?diag| ~ 0.2 for 
L > 6. With increasing L, this plateau gradually gets higher but no side peaks show up. Since the normalization 
condition of P s (<?diag) inhibits the plateau of finite width getting higher indefinitely, one plausible asymptotic behavior 
of -P s (?diag) might be that it eventually converges in the L — > oo limit to the Gaussian distribution around Qdiag = 0. 
In that sense, the observed behavior is consistent with the spin disorder at this temperature. However, solely from 
the present data, we cannot completely rule out the possibility that the ±qf A /3 peaks characteristic of the SG LRO 
eventually show up for still larger L. Anyway, in the range of sizes studied L < 10, we have observed no sign of such 
side peaks, in contrast to the 5D case where the side peaks appear already for L = 5. Another possibility may be 
that the 4D model exhibits a finite-temperature SG transition but the SG state is a critical phase with a vanishing 
SG order parameter qf A = 0, as expected for the system at its LCD. However, as argued below, such a LCD behavior 
is not supported from our data of the critical properties. 
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FIG. 6: Overlap distribution function of the diagonal part of the spin, in 4D at T/J = 0.26 (a), in 5D at T/J — 0.40 (b), and 
in the SK model at T/J = 0.10 (c). 



F. Chiral-glass susceptibility 



From the data presented above, we have concluded that in 4D the chirality is an order parameter of the transition, 
but not so in 5D and in the SK model: In the latter cases, the order parameter of the transition is the spin, while 
the chirality order is parasitic to the spin order. In order to examine further the validity of such a picture, we show 
in Figs0 the temperature and size dependence of the reduced chiral-glass susceptibility, where \ x is plotted as a 
function of the reduced temperature (T — Tcg)/Tcg on a log-log scale. As determined above, Tqg is taken to be 
T CG /J = 0.38 in 4D and T CG /J = 0.60 in 5D. In the SK case, we put T CG /J = \ which is exact. 

As can be seen from Figs[7| (b) and (c), the reduced chiral-glass susceptibilities of the 5D and the SK models do 
not exhibit any sign of divergence: They stay small at any temperature T > Tcgj and most notably, \x S e ^ s smaller 
with increasing the system size L. Such a size dependence is completely opposite to the one expected for a diverging 
quantity in the critical region. We note that even at temperatures close to the transition temperature, no sign of the 
reversal of the size dependence is discernible. In fact, in the SK case, such a non-diverging behavior of x x is just as 
expected. In the SK model, the standard SG exponents are known to be a = — 1, /?sg = 1 and 7sg = 1- Then the 
chiral-glass exponent /3cg is expected to be /3cg = 3, because the chirality is cubic in the spin variables. Then, the 
chiral-glass susceptibility exponent is obtained as 7cg = — 3 from the scaling relation a + 2/?cg + 7cg — 2. Negative 
7cg means that the chiral-glass susceptibility of the SK model does not diverge at the SG transition. Very much 
similar behavior observed in 5D suggests that the the chiral-glass susceptibility of the 5D model does not diverge at 
the transition, either. Hence, our observations for x x are fully consistent with our previous finding that the order 
parameter of the transition in 5D and in the SK model is the spin, not the chirality. 
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By contrast, in 4D, \ x exhibits a different behavior. Although in the investigated temperature regime Xx stays 
rather small and tends to decrease with increasing L, similarly to the behavior observed in 5D and in the SK limit, its 
size dependence is about to change in a close vicinity of Tqg- More specifically, the L = 10 data catch up the L = 8 
data at t ~ 0.04, and at temperatures further close to Tec exceeds the L = 8 data, where the data show significant 
finite-size rounding preventing the observation of the asymptotic critical behavior. This suggests that the critical 
region of the chiral-glass transition might be very narrow in 4D, limited to the regime t 5 10~ 2 . In the temperature 
range outside this, the chirality exhibits a mean-field-like non-diverging behavior similar to the one of the SK model. 
Although we cannot directly get into this narrow critical region in our present simulation due to the computational 
limitation, the observed behavior of \x °f the 4D model indeed suggests that such a crossover from the mcan-field-likc 
behavior at t ^ 10~ 2 to the diverging critical behavior at t 1$ 10~ 2 does occur. Unfortunately, inaccessibility to the 
asymptotic critical region prevents us from estimating the chirality exponents. In order to estimate the asymptotic 
chiral critical exponents, one needs to approach the temperature regime t IS 10 -2 with larger lattices L > 10, which 
is not feasible with the computational capability presently available to us. 
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FIG. 7: Temperature and size dependence of the reduced chiral-glass susceptibility in 4D (a), in 5D (b), and in the SK model 
(c). 



G. Critical Exponents 

In this subsection, we analyze the critical properties of the SG transition in each case of 5D and of the SK model. 
In the case of 4D, our analysis has suggested that the transition is a pure chiral-glass transition, the associated chiral 
critical regime being narrow, t IS 10~ 2 , which prevents us from estimating the chiral critical exponents. 
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In 5D, we have concluded that the transition is the standard SG transition, i.e., the order parameter is the spin, 
not the chirality, although the chirality also takes a nonzero value in the SG ordered state reflecting the noncoplanar 
character of the spin order. We estimate the associated SG exponents via the standard finite-size-scaling analysis of 



.,(2) 



based on the relation, 



(2) 



L 



-{d— 2+J7sg) ffJ^/vsG 



\T-Tt 



SG 



(21) 



where the T§g value is set to the best value determined above, Tgc/ J — 0.60. The best estimates of the SG exponents 
are z/go — 0.6(2) and ?7sg = —0.8(2): See Fig|SJ From the standard scaling relation, we get other exponents as 
a = —1.0(5), ftsG = 0.7(3) and 7sg = 1-7(5). One sees that these exponents are not far from the mean-field exponents 
expected above the upper critical dimension D = 6. 

In the SK case, mean-field exponents should be exact. Indeed, as shown in FigEl our data of qi 2 ^ are entirely 

25.26 



consistent with such a mean-field behavior, (3sg = 1 and ?7sg 

In concluding this section, we touch upon the near-critical behavior of the spin observed around the chiral-glass 
transition point in 4D. Although our data of the Binder ratios and the overlap distribution functions given above have 

(2) 

strongly suggested that the Heisenberg spin remains paramagnetic even below Tqg in 4D, the present qs data can be 
scaled reasonably well with assuming T^q/ J — Tqg/J — 0.38(2), at least from the purely numerical viewpoint. Such 
a constrained finite-size scaling analysis of qs yields the estimates i>' SG = 1-3(2) and rj' SG = —0.7(2): See FiellOl (a). 

Note that the value of v is far from the LCD value, v = oo. We also note that the same q^ data can also be fitted 
with a comparable quality by assuming a zero-temperature SG transition, Tg G = 0: See Fie llOl (b). As mentioned, 
we believe that the is' SG and rj' SG values obtained by assuming T$q — T G g do not represent true asymptotic exponent 
values, but just represent effective exponents describing the short-scale near-critical phenomena which is an echo of 
the chiral-glass transition. Indeed, at short scales, the chirality is never independent of the spin by its definition. 
Hence, the behavior of the spin-correlation related quantities might well reflect the critical singularity associated with 
the chirality up to certain length and time scales. 
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FIG. 8: Finite-size scaling plot of the spin-glass order parameter in 5D. 



V. SUMMARY AND DISCUSSION 



In summary, we performed a large-scale equilibrium MC simulation of the 4D, 5D and SK Heisenberg spin glasses. 
In 5D, the model exhibits a single SG transition at a finite temperature, reminiscent to the one of the corresponding 
mean-field model. Below the transition temperature T^q/J — 0.60(2), the spin is frozen in a spatially random 
noncoplanar configuration. Although the SG order accompanies a finite chiral-glass LRO reflecting the noncoplanar 
nature of the spin order, the order parameter of the transition is the spin, not the chirality, and the chiral-glass 
susceptibility remains nondiverging at T — Tsg- Similar behavior is also observed in the Heisenberg SK model. 
The SG exponents in 5D are estimated as a = —1.0(5), /?sg = 0.7(3), 7sg = 1-7(5) and ^sg = 0.6(2), most of 
which are rather close to the mean-field exponents. Since the upper critical dimension of the SG is believed to be 
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six, the observed closeness to the mean-field values seems reasonable. Our data suggest that the SG ordered state 
accompanies a peculiar phase-space structure, namely, a one-step-like RSB, at least in its chiral sector. Such a one- 
step-like character of the RSB is at variance with the full (hierarchical) RSB realized in the SK model corresponding 
to D = co. It means that the RSB pattern of the ordered state changes its nature at some borderline dimensionality, 
presumably at the upper critical dimension D = 6. 

In 4D, the model exhibits a significantly different behavior from the 5D and the SK models. Bulk of our data, 
particularly the Binder ratio, indicate that the 4D model exhibits a pure chiral-glass transition at a finite temperature 
Tcg = 0.38(2) without accompanying the standard SG order. At the chiral-glass transition, however, the spin becomes 
almost critical which manifests itself as a pseudo-critical phenomenon observable at short length scales. The critical 
region associated with the chiral-glass transition is very narrow, limited to t ^ 10~ 2 , suggesting that the dimension 
four is close to the marginal dimensionality. The SG transition occurs either at Tsg — or at a finite temperature, 
but below the chiral-glass transition temperature, Tsg < Tcg- Our data suggest that the chiral-glass ordered state 
accompanies a one-step-like RSB in the chiral sector. Again, such a one-step-like character of the RSB differs from 
the full (hierarchical) RSB realized in the SK model. 

Next, we wish to compare out present results on the 4D and 5D Heisenberg SGs with those of the previous authors. 
To our knowledge, our results for the chiral order are new. Concerning the spin order, our present conclusion, i.e., 
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the presence of the SG LRO in 5D and the absence of it in 4D, is consistent with most of numerical simulations, in 
particular, with the one of Stauffcr and Binder—. Our conclusion, however, is at variance with that of Coluzzi, who 
suggested that the SG LRO set in at a finite temperature T$g/J — 0.5^. Although the numerical data themselves 
seem to be consistent between the two works, Coluzzi simulated rather small lattices L < 5 and high temperatures 
T/J > 0.5, which hampered a direct examination of the asymptotic ordering behavior. Our new data for larger 
lattices L < 10 and for temperatures including lower ones, T/J > 0.26, have clarified that the transition occurs in the 
chiral sector at Tqg/J — 0.4, which is somewhat lower than Tsg estimated in Rcf. 20 . Furthermore, the Heisenberg 
spin appears to remain paramagnetic at the chiral-glass transition point on sufficiently long length scales, i.e., the 
transition at Tcg/J — 0.4 is not the conventional SG transition, but a pure chiral-glass transition. 

Finally, we wish to discuss implication of our present result to the 3D case. The behavior of the 4D model observed 
in the present work is qualitatively similar to the one of the 3D model observed in Refsii2*i^, except that the behavior 
of the 4D model looks much more marginal. For example, the reduced chiral-glass susceptibility of the 3D model is 
much larger in magnitude than that of the 4D model, and the associated chiral-glass critical region is much wider in 
3D than in 4D. As one judges from the size dependence of the reduced chiral-glass susceptibility shown in Fig. 3 of 
Ref4£, the width of the chiral critical region is about 10 _1 , which should be compared with our present estimate for 
the 4D model, 10 -2 . All these suggest that the spin-chirality decoupling is more eminent in lower dimensions. As the 
dimensionality is increased, the spin-chirality decoupling tends to be suppressed. In 4D, the spin-chirality decoupling 
still seems to persist, but it is limited only to a very narrow temperature region close to the transition temperature, 
suggesting that 4D is close to the borderline dimensionality. As the dimensionality is further increased, the spin- 
chirality decoupling no longer occurs. There, the order parameter of transition is the spin, not the chirality. This 
is indeed the case for 5D. However, at least in the case of 5D, the SG ordered state exhibits a peculiar one-step-like 
RSB, which differs in character from the full RSB of the D — oo SK model. Estimated SG critical exponents of the 
5D model are rather close to the mean-field values, which is consistent with a common belief that the mean-field SG 
exponents arise above six dimensions. 



The numerical calculation was performed on the HITACHI SR8000 at the supercomputer system, ISSP, University 
of Tokyo, and Intel Pcntium4 1.8GHz PCs in our laboratory. The authors are thankful to Dr. K. Hukushima and Dr. 
H. Yoshino for useful discussion. 



In this appendix, we give the derivation of Eq.Q. It describes the self-overlap part of the diagonal spin-overlap 
distribution function in the thermodynamic limit, when the SG ordered state with a nonzero EA order parameter 
qf A > exists. Since the diagonal spin-overlap <?diag transforms nontrivially under global 0(3) rotations, even a 
self-overlap part of the distribution function is not just a simple delta function located at <?diag = ±9? A , but exhibits 
a nontrivial behavior given by Eq.J7|). 

We consider a diagonal spin-overlap between a particular spin state described by the configuration Si and a state 
generated from this state via a global 0(3) rotation R, 



We first consider the case of proper rotations with det(iZ)=l. The 50(3) rotation matrix R is known to be 
parametrized by the Euler angles, O and as 
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APPENDIX A: DERIVATION OF EQ.© 




(Al) 




cos O cos $ cos ^ — sin $ sin \I/ 
— cos O cos $ sin ^ — sin <I> cos ^ 
sin O cos $ 



cos O sin $ cos ^ + sin $ cos — sin O cos \ 
— cos O sin $ sin + cos $ cos ^ sin O sin 
sin O sin $ cos O 




Then, (jviiag can be written as 
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The spin direction at each site can be represented as Si = (sin 8i cos fa, sin 6^ sin fa, cos 4>i). If one notes the fact that, 
in the SG ordered state, the spin direction is entirely random on long length scales giving a uniform distribution on 
a sphere in spin space, one can replace in the thermodynamic limit the summation over spins by the integral over 
spin directions as (1/^0X^1 ~~ * (V^ 71 ") Jj^dcosfl d<p. Then, only the diagonal terms containing survive in 
Eci. (|A2(> . leading to 

cos(3> + *) + 1 ^ cos($ + vp) cos($ + vp) + 1 cos($ + <£) . , , 

9dia g = - ^ cos6 + '- = * ^ x + K — '- , (A3) 

where x = cos 0. Note that this is a function of the rotation matrix R only, not depending on the spin configuration Si 
any more. The overlap (jdiag takes various values depending on the 0(3) matrix R. We then consider the distribution 
of 9diag arising from the distribution of R, or equivalently, x, $ and The appropriate measure is —1 < x < 1, 
< $ < 2ir and < \t < 2ir being uniform. It is convenient to change the variables from ($, fy) to (a, (3) = 
(<f> + (-$ + *)/2), where < a < An and < (3 < n. With this change of the variables, Ea . (| A3|) becomes 
independent of /3, and is given by 

cos a + 1 cos a , . . 

9diag = g x + — — . (A4) 



The distribution function -P s (<Zdiag) is proportional to 

f S («diag) oc f dx da 
J dg d iag 



d(?diag 

Qic(gdlag) 3 



1 + cos a 



da , (A5) 



where we have used Ea. (|A3(l . Note that, for a given gdiag, the integral with respect to a is restricted to the range 
[0,a c ], with a c (gdiag) = cos _1 [(3qdiag — l)/2]- This may be seen from Fig llll where we plot gdiag as a function of x 
for various a. Obviously, for a given gdiagj n o contribution to the integral arises from the region of a between [a c , n]. 
The Qdiag-dependence of -P s (<Zdiag) arises from this upper limit of the integral. The integration in Eci. (|A5(l can be easily 



9dia g (a;,a!) 




-1/3 



FIG. 11: Sketch of (/diag as a function of x for various a. For a given value of (faiag, there is no possible value of x for a range 
of a satisfying a c (gdia g ) < a < n. 



carried out to yield, 

-P s (<7diag) oc tan 



«e(ffdiag) _ „l/ 2 / 1 ~ gdia 



3(/diag + 1 



So far, we have considered proper rotations. The contribution from improper rotations, which can be viewed as 
proper rotations combined with the spin inversion Si — > —Si, may be obtained immediately by the replacement (faiag — * 
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— 9diag- Adding the contributions from both proper and improper rotations with equal weights, and reproducing the 
appropriate normalization factor, we get 

Finally, we note that at finite temperatures a state should be regarded as a pure state. The spin length is then no 
longer unity, and the unity in Ea. ljA6|) should be replaced by qf A . We finally obtain 

3 toag + qf A + V -3, diag + qf* J • ( A? ) 

which is Eq.Q. 

The derivation above, valid in the thermodynamic limit N — ► oo, is quite general. In order to get some feeling 
about the finite-size effect, we also compute -Ps('Zdiag) numerically for finite- N Hcisenberg spins, the direction of which 
is assumed to be completely random in three-component spin space. More specifically, we prepare a random and 
uncorrelated configuration of N spins, numerically generate 0(3) rotation matrix R with appropriate measure (i.e., 
the one generated randomly from the uniform distribution of— 1<x<1,0<<1 > < 2tt, < \& < 2ir and the 
determinant ±1), operate R to the initial spin configuration, and compute the diagonal spin-overlap (faiag between 
the initial and the 0(3)-rotated spin configurations. We generate 10 4 distinct 0(3) matrices for a given initial spin 
configuration, and generate several hundreds of initial spin configurations, -P s (<Zdiag) being accumulated over these 
procedures. The result is shown in Fig^| The N = oo result analytically obtained above is also shown. It can be 
seen that the rounded peak at ± ^qf A grows as N increases, eventually exhibiting a divergent behavior in the N = oo 
limit. Of course, the finite-iV result computed here is valid only for non-interacting spins. It would differ from the 
corresponding result for the interacting system, in contrast to the analytical N = oo result which is valid even for the 
interacting system. It still gives some feeling how the finite-size rounding takes place in finite- N SG models. 
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